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Abstract 

Shape Dynamics is a formulation of General Relativity where refoliation invariance is traded for local 
f— jh spatial conformal invariance. In this paper we explicitly construct Shape Dynamics for a torus universe 

in 2+1 dimensions through a linking gauge theory that ensures dynamical equivalence with General 
Relativity. The Hamiltonian we obtain is formally a reduced phase space Hamiltonian. The construction 
of the Shape Dynamics Hamiltonian on higher genus surfaces is not explicitly possible, but we give an 
explicit expansion of the Shape Dynamics Hamiltonian for large CMC volume. The fact that all local 
constraints are linear in momenta allows us to quantize these explicitly, and the quantization problem for 
Shape Dynamics turns out to be equivalent to reduced phase space quantization. We consider the large 
CMC-volume asymptotics of conformal transformations of the wave function. We then use the similarity 
of Shape Dynamics on the 2-torus with the explicitly constructible strong gravity (BKL) Shape Dynamics 
Hamiltonian in higher dimensions to suggest a quantization strategy for Shape Dynamics. 



1 Introduction 



Shape Dynamics is constructed as reformulation of General Relativity in which spacetimc refoliation in- 
variance is traded for spatial conformal invariance that preserves the total spatial volume [1, 2, 3] using 
of a linking gauge theory. The development of Shape Dynamics was inspired by Dirac's work [4] on CMC 
(constant mean extrinsic curvature) gauge, York's method for solving the initial value problem [5, 6, 7] and 
Machian ideas developed by Barbour and collaborators [8, 9]. The explicit construction of the Shape Dynam- 
ics Hamiltonian however requires the general solution of a partial differential equation, which is equivalent 
to partially solving the initial value problem of General Relativity using York's method. This complication 
obstructs many straightforward investigations. To learn about Shape Dynamics it is therefore very valuable 
to consider exactly solvable nontrivial gravitational models. This provides the main motivation for this 
paper: We consider a nontrivial model in which Shape Dynamics can be constructed explicitly allowing us 
to study its generic features. 

Probably the best known example of a nontrivial exactly solvable gravitational system is pure gravity 
on the torus in 2+1 dimensions [10, 11, 13, 12]. The technical reason for the simplifications in this model is 
two-fold: First, one is able to solve the initial value problem of ADM gravity explicitly on the 2-torus. This 
is important for the construction of classical Shape Dynamics and occurs only on the 2-torus and 2-sphere; 1 
pure gravity on higher genus surfaces is more intricate since we lack methods to solve the initial value problem 
in general. Second, the reduced phase space (after solving for initial data) is finite dimensional, which is 
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1 The sphere is a degenerate case, since it admits only one canonical pair of degrees of freedom (the volume and the mean 
extrinsic curvature). This admits only the de Sitter solution which contains no interesting dynamics. 



a generic feature of pure gravity in 2+1 dimensions. This is important for quantization, because a finite 
dimensional system admits generic quantum theories while nontrivial quantum systems with infinitely many 
degrees of freedom are sparse. 

The plan for this paper is as follows: We start with the explicit construction of pure Shape Dynamics 
on the 2-torus in section 2 and show its equivalence with General Relativity using the method of linking 
gauge theories. The trading of rcfoliation invariance for local spatial conformal invariance turns all local 
constraints into phase space functions that are linear in the momenta, while the remaining Shape Dynamics 
Hamiltonian turns out to formally coincide with the reduced phase space Hamiltonian, which at large CMC- 
volume becomes the conformal constraint that changes the total volume. In section 3 we consider Shape 
Dynamics on a higher genus surface, which does not admit straightforward solutions to the initial value 
problem. We attack this problem by constructing a perturbation expansion and we recover a fully conformal 
theory in the large volume regime. We find that the generically nonlocal Hamiltonian becomes the integral 
over a local density at large volume and turns again into the conformal constraint that changes the total 
volume. We then use the classical results to quantize Shape Dynamics on the torus in section 4. Due 
to linearity of the local constraints, one can implement them at the quantum level and thus quantize the 
analogue of the Wheeler-DeWitt equation, i.e. infinitely many degrees of freedom. In an excursion we 
discuss the large-CMC-volume asymptotics of Shape Dynamics wave functions, which resembles aspects of 
a quantum gravity/CFT correspondence. Before we conclude in section 6 we discuss which results of this 
paper could generalize to higher dimensions and which are specific to 2+1 dimensions in section 5. 

2 Equivalence of General Relativity and Shape Dynamics 

In this section we establish the equivalence between General Relativity and Shape Dynamics on the 2+1 
dimensional torus universe by explicitly constructing the linking theory relating the two. For simplicity we 
assume a positive cosmological constant A. We start with the general construction of the linking theory 
before focusing on the torus, which allows us, other than in higher dimensions or even on a higher genus 
surface, to explicitly work out Shape Dynamics. 

A convenient way to construct the linking theory is to follow the "best matching procedure" [1] starting 
with the ADM Hamiltonian on the usual ADM phase space Tadm, expressed in terms of the metric g a b(x) 
and its canonically conjugate momentum density w ab (x): 

H = S(N) + H(0, 
S(N) = Jd 2 xN(-L= i TT ab G abcd n cd -^\i\(R-2A)Y (1) 

H(0 = fd 2 xir ab £ igab . 

Here G a b c d = 9ac9bd — 9ab9cd is the 2-dimensional super-metric, and S and H denote the ADM scalar and 
diffcomorphism constraints. To best match with respect to conformal transformations that preserve the total 
volume we consider the ADM phase space as a subspace of a larger phase space T ext = Tadm x F^, where 
r^ is the phase space of a scalar field <fi(x), whose canonically conjugate momentum density is denoted by 
7T0(x). The phase space functions on Tadm are naturally identified with those phase space functions on T ext 
that arc independent of (<fi, ir^). We can thus recover usual ADM gravity in this larger system by introducing 
an additional first class constraint 

Q{x) := 7r (z) « (2) 

and add it smeared with a Lagrange multiplier p(x) to the ADM Hamiltonian, which is now H = S(N) + 
H(£) + Q{p). Let us now consider a canonical transformation from (g a b,^ ab , (f), tt^) to (G a b,TL ab ,^,Tl ( f,) 
generated by the generating functional 



F= d 



2 x (g ab {x)e 2 ^T\ ab {x) + 0(x)II^(x)) • (3) 

Here <p is defined in terms of <fi by subtracting a spatial average, that has a non-trivial dependence on the 
metric, 

fa):=<j>{x)-\\*{<?*). (4) 



where we use the shorthands (f) g = V 1 J d 2 xy/\g\f and V g = J d 2 xy/\g\. Notice that we constructed <f> such 

that the conformal factor e 2 ^ preserves the total volume. The canonical transformation of the elementary 
variables can be worked out explicitly: 

G ab (x) = e 2 ^g ab (x), 

U ab (x) = e~ 2 ^) (Tr ab (x) - ±T/\i\(x)g ab (x)(n) (l - e 2 ^)) , 

*(*) = 0(a), 

II^(x) = ■Kj,{x)-2 (k(x) -(Tr)y/\g\(x)J , 

using shorthand notation 7r(x) = ir ab (x)g ab (x) and (7r) = V -1 J d 2 xir(x). This transformation leads us to 
the constraints of the linking theory: 

H = S(N) + H(0 + Q( P ), 
S(N) = fd 2 xN[^(w ab G abcd n cd -l(w-{7T)(l-e^)^j 2 + l7r 2 

-7l5|(%]-2A^-2Ae 2 ^)], (6) 

H(0 = fd 2 xe- 2 ^ U ab (x) - \^\g~\{x)g ab {x){n) (l-e 2 ^)) (c s e 2 ^g\ (x), 

Q(p) = fd 2 xp(x)(TT4x)-2(n(x)-(n)^\i\(x))). 

One can check that after integrating by parts and using Q — the constraint H(£) turns into the usual form 
of the diffcomorphism constraint H(£) — J d 2 x {Tr ab C^g a b + ^<pC-^<p) ■ We will use this form of the constraint 
below. The linking theory thus contains the usual diffcomorphism constraint, a conformal constraint that 
preserves the total 2-volumc and a scalar constraint that arises as a modification of the ADM refoliation 
constraint. 

2.1 Linking Theory on the Torus 

We will now exploit some special properties of two dimensional metrics on the torus to simplify the constraints 
(6). First of all, it is a well-known fact that all metrics on the torus are conformally flat. The space of flat 
metrics modulo diffcomorphisms is finite dimensional and admits a convenient paramctrization of the space 
of metrics on the torus. We follow [12, 13] where possible. 

To make this more explicit let us fix a global chart on the 2-torus T 2 , which allows us to uniquely identify 
any point x G T 2 with its coordinates (x 1 ,^ 2 ) £ [0, l) 2 . In these coordinates we can decompose an arbitrary 
metric g ab as 

g ab (x) = e 2X ^(f-g) ab (x), (7) 

where A is a conformal factor, / is a (small) diffcomorphism and g a flat reference metric. We can make this 
decomposition unique by requiring g to be of the form 

g = — (dx 1 <g) dx 1 + r x {dx l <g> dx 2 + dx 2 ® dx l ) + (r 2 + r%)dx 2 <g> dx 2 ) (8) 

where r = (ti, ti) denote the Teichmiillcr parameters. There is a slight redundancy left in the decomposition 
having to do with the fact that / is only determined up to an isometry of g, i.e. up to translations in x 1 and 
x . If we require / to leave (0, 0) invariant, we really obtain a one to one map between metrics on the torus 
and the data (A, /, r). 

Wc can explicitly decompose the momentum density 2 

Tr ab = e- 2X fp ab + -irg ab + y/\§\ (D a Y b + D b Y a - g ab D c Y c ) j , (9) 



"All indices here are raised with the reference metric g 



in terms of a trace 7r, a vector field Y and a transverse traceless tensor density (w.r.t. g), which we can 
explicitly parametrize by 

P = o ((( r i 2 + r 2 )P 2 ~ 2r i T 2Pi) 9i ® (9i + (r 2 pi - np 2 ) (e?i <g> <9 2 + 9 2 (g> 3i) +p 2 di % d 2 ) ■ (10) 

This decomposition is such that 7r is conjugate to 2A and pi is conjugate to r^. 
Writing the linking theory constraints (6) in terms of A, /, r, 7T, V and p we get 



S(iV) = Jd 2 xN 



(V b + (PY ab )^/\g~\)g ac g bd (p« l + (PY cd )^\f\) 



(11) 



-i(7T - (tt) VM(1 - e 2 ^+ A ))) 2 ) - VM (fc - 2A(0 + A) - 2Ae 2 (^+ A )) 
H(0 = /d 2 ^ a (^|Ar a + ie 2A 5 a ( e - 2 V)+^0 ja ) 
Q(p) = fd2 X p(x)^(x)-2(n(x)-(7r)yW\(xy^)), 

where we used the shorthand PY a {, = DaXb + DbY a ~ gabg cd D c Yd- 

To complete the definition of a linking theory, wc specify two sets of gauge-fixing conditions, 

<j>(x)=0 for GR and n^x) = for SD, (12) 

which wc will now use to reconstruct General Relativity and Shape Dynamics respectively. 

2.2 Recovering General Relativity 

To recover standard ADM gravity on the torus let us impose the gauge-fixing condition <j){x) = to the 
linking theory. To perform the phase space reduction from the extended phase space to the ADM phase 
space, we need to fix the Lagrange multipliers such that the gauge-fixing is propagated. Since the momentum 
density 7r^ occurs only in the constraints Q, we have to solve 

= {Q(p),<t>(x)}=p(x) (13) 

for the Lagrange-multiplier p, which implies p = 0, hence the constraints Q(x) are gauge-fixed and in fact 
drop out of the Hamiltonian, which becomes independent of tt^. Hence, one can perform the phase space 
reduction by setting <f) = 0, p = and 7T0 arbitrary 3 in equation (6). The Hamiltonian on the ADM phase 
space thus reads 

H = S(N) + H(0, (14) 

where S(N) and H(£) are precisely the scalar- resp. diffeomorphism constraint of General Relativity in the 
ADM formulation as given in equation (1). We note that we explicitly retained refoliation invariance. 

2.3 Recovering Shape Dynamics 

To recover Shape Dynamics we employ the gauge- fixing condition itj,(x) = 0. We will see that the decomposed 
form (11) of the constraints allows us to find the explicit Shape Dynamics Hamiltonian through a phase space 
reduction (0, 7r^) — > (</>o, 0). To find this map we can use tt^ = 0, so the Q constraints become 

Q(p)= J d 2 xp(x) (n(x) - (tt) y/Mx)) , (15) 

which implies that w(x) is a covariant constant. Using this and 7T0 = 0, we find that the diffeomorphism 
constraint implies that 

PY ab (x) = 0, (16) 



J Had the constraints Q(p) not dropped out after gauge fixing p, wc would have had to solve Q = for tt^ to complete the 
phase space reduction. 



which implies that the scalar constraint is independent of Y a (x). Using these simplifications in the scalar 
constraints, we find 

S(N) = J d 2 x^\§\N f e -2(i+x) ~9ac9bdP^ b P cd _ l e 2(0+A) (W 2 _ 4 A) + 2A(0 + A)) . (17) 

Using (8) and (10) we find that 

g a c9bdP ab p cd r| 2 2 

j^j = y(Pi+P 2 ) (18) 

is a spatial constant in the chosen chart. 



We see that the constraints S(N) would be solved if we were able to choose e 4( -* +A - ) - 



|g|«7r)2-4A) ■ 

However, this is in general obstructed by the volume-preservation condition J d 2 x^J\§\e 2 ^ +x > — V . This 
means that the constraints generating the rcfoliations arc not completely gauge fixed by the condition 
tt^,(x) — 0. Indeed it turns out that among the infinitely many constraints S(N) one remains first class, 
which after phase space reduction becomes our Shape Dynamics Hamiltonian Hsd- More concretely, there 
exists a lapse Nq such that S(Nq) Poisson commutes with 71^, i.e. it satisfies the lapse fixing equation 



{5(iVo),7r (x)} = F No (x)-e^+^^/\g~\(F No )=0 

where F N = N (-2e~ 2 ^+ x ">g ac g bd p ah p cd - \g\e 2 ^+ x \{n) 2 - 4A)) + ^/\f\AN. 



If one imposes on No a normalization condition J d 2 x\J\§\e 2 ^ +x ^ Nq = V, then (19) has a unique solution. 
We want to project out this first class part S(N ) from the full set of constraints S(x) to end up with a 
purely second class set of constraints S(x) that we can solve. We can perform the projection in different 
ways, but a particularly convenient way of doing this is by defining 

S(x) := S{x) ^M/^)e 2 ^> + ^», (20) 

which automatically satisfies S(No) = 0. Identifying Hsd = S(Nq), we arrive at the modified Lichnerowicz- 
York equations 

= S(x) = VM(e- 2( ^ +A) g °^f| aVd ~ 2-eW>(W 2 -4A + 2«p) - 2A(^ + A)) , 
V = J d 2 xy/\f\e 2 ^+ x \ 

which we need to solve for <\> and Hsd- A solution is found by taking <f> + A to be spatially constant. More 
precisely, from the second equation it follows that 

4> = -X + -lnV. (22) 

Now Hsd can be easily determined from the first equation in (21), 

1 7t ri -n a b^cd T/ ^-2 T/ 

Hsd = l gncgt g] P - \ (W 2 - 4A) = ^{p\ + p 2 2 ) - \ ((ir) 2 - 4A) . (23) 

Notice that to find Hsd = S(Nq) we did not have to solve the lapse fixing equation explicitly. In this case 
we can solve (19) straightforwardly using the fact that <f> + A is constant and the result is simply No = 1. In 
general however the lapse fixing equation is quite complicated and we are lucky that we don't actually have 
to solve it to derive the Shape Dynamics Hamiltonian (as we will again see in section 3). As a matter of 
fact, as a constraint S = is completely equivalent to 



S(x) - (S) y/\§\(x)e 2 ^^ +x ^ = 0, (24) 

which does not refer to a lapse at all. 



The Shape Dynamics Hamiltonian Hsd (23) is exactly the reduced phase space Hamiltonian constraint. 
The more familiar Hamiltonian -ffyork generating evolution in York time (ir) (see e.g. [13] section 3.3) is 
obtained by noting that the variable canonically conjugate to (it) is V and therefore by solving Hsd = 0, 



ff York = F = r 2 -pp2L (25) 

We can now perform explicitly the phase space reduction of the linking theory and describe Shape 
Dynamics on the ADM phase space through its total Hamiltonian and first class constraints 

H = JVH SD + H(0 + C{p) 

Hsd = 4(p 2 + P 2 )-^((ir) 2 ~4A) 

H(0 = Jd 2 xir ab ^g ab _ ( 26 ) 

C{ P ) = f<Pxp(*-{*),/\g\)- 

The gauge symmetries are indeed spatial diffeomorphisms, conformal transformations that preserve the total 
volume and global time reparametrizations. Despite the different set of symmetries, the equivalence with 
standard General Relativity is obvious: The Shape Dynamics Hamiltonian coincides on the reduced phase 
phase with the CMC Hamiltonian, while the constraints C provide the CMC gauge-fixing conditions. 

Although we know the Shape Dynamics Hamiltonian explicitly on the torus, it is instructive to observe 
that the Shape Dynamics Hamiltonian constraint H can be expanded in powers of the inverse volume, 
because it shows two properties that we can investigate in more complicated models. This expansion is a 
systematic approximation to Shape Dynamics that is a good approximation in an asymptotic large volume 
regime, i.e. where V — > oo while keeping the other degrees of freedom finite. In this regime we find two 
important features of Shape Dynamics: 

1. Asymptotic Locality: The leading order of the Hamiltonian, which becomes exact in the limit 
V —> oo, is (ir) 2 — 4A + 0(V~ 2 ) « 0. As a constraint, this is equivalent to 

V ({it) - 2\/A~) = J d 2 x (n(x) - 2VAy/\g~\(xj) « for V ->■ oo (27) 

which is diffcomorphism invariant as the integral over a local density and by inspection invariant under 
conformal transformations that preserve the total volume. 

2. Full Conformal Invariance: Since the Shape Dynamics Hamiltonian constraint is asymptotically 
equivalent to (it) — const, w 0, we can add it to the conformal constraints C to obtain in the large 
volume limit C(x) + Hsd = tt(x) — const, w 0, which generates full conformal transformations, i.e., 
including those that change the total spatial volume. Notice that this requires us to interpret the 
Shape Dynamics Hamiltonian as a constraint, rather than a generator of physical dynamics. 

Let us have a quick look at the 2-sphere: The linking theory and phase space reduction can be performed 
following the same steps as on the torus with two small modifications. 1) there are no Teichmiillcr parameters 
on the sphere, so there is only one canonical pair of physical degrees of freedom (V and (it)); 2) the total 
spatial curvature does not vanish, but is 87r. Shape Dynamics on the sphere thus takes the form of equation 
(26), except for the Hamiltonian, which is Hsd = — y(( 71 ") 2 — 4A) — 87r. 

3 Higher Genus Surfaces 

On a higher genus surface, we can still use the decomposition analogous to equations (7) and (9), but the 
explicit construction of the Shape Dynamics Hamiltonian constraint on the torus rested on the explicit 
solvability of the modified Lichnerowicz-York equation (21). The Lichncrowicz-York equation on a higher 
genus surface (or in higher dimensions) is not explicitly solvable. We thus restrict our construction of Shape 
Dynamics to an approximation scheme and consider an expansion that becomes exact in the large volume 
limit. 4 Again, we follow [12, 13] when possible. 



J Scc [14] for an analogous expansion. 



3.1 Preparations 

In genus g > 2 we can perform a decomposition analogous to (7) and (9), 



g ab (x) = e 2 ^(f*g) ab (x), 

TT ab (x) = e- 2X [p ab (x)+ 1 1 g a \x)Tr(x) + ^/W\(x)g ab (x)g cd (x)PY cd (x)), (28) 



where now we take the reference metric g to be of unit volume / d 2 X\J\g\ = 1 and constant scalar curvature 
R. According to the Gauss-Bonnet theorem, R is given by 

R = -8n(g-l). (29) 

Modulo diffcomorphisms the space of such metrics corresponds to the genus g Teichmuller space, which has 
dimension 6g — 6. Unfortunately no simple explicit parametrization for g is known, so we will keep the 
parametrization implicit. 

We can again write the linking theory using the decomposition (28). The only difference compared to 
the constraints (11) for the torus is the subtraction from S(N) of a spatial curvature term. 

When wc impose the gauge fixing tt^ = 0, we obtain the analogue of (17): 

S(N)= f d 2 x^/\f\N (V 2 ^+ A > ^9bdP^P cd _ 1 e 2(0+A) ((?r) 2 _ 4A) +2 A(4> + \)-rY (30) 

Reusing our discussion for the torus, we construct the second class part S of S according to (24), 



S(x) = S(x) - (S),/\§\(x}e 2 M a ') +x (*». (31) 

Identifying the remaining first class constraint (S) with Hsd/V, we obtain the modified Lichnerowicz-York 
equations for genus g > 2, 



= S(x) = y/\f\ ( e ~2(0+A) g^|!Vf _ l e 2(0 + A) ( ^ )2 _ 4A + 2 flp } + 2A( ^ + X)-R), 

V = J^xy/Wle 2 ^^- 



To simplify the notation let us define /i = (f> + X — ^hiV and p 2 := 9ac9bd P — E — . Then (32) can be written 

^p 2 e- 2 ^±v((7r) 2 -4A+^^)e 2 » + 2Av-R = and (e 2 »)- g = 1. (33) 

In the following wc will drop the subscript g and keep in mind that averages (•) are taken with respect to g 
(except for (it)). 

Equation (33) is nearly identical to the standard Lichnerowicz-York equation in 2+1 dimensions. The 
only difference is that we have a restriction on /i and to compensate this we have an additional constant 
Hsd to solve for. The existence of a unique solution for /i and Hsd (as a function of g a i,, p ab , V and (tt)) 
is a direct consequence of the existence and uniqueness properties of the usual Lichnerowicz-York equation. 

The key simplification that allowed us to explicitly construct Shape Dynamics on the torus is that there 
one can choose the constant curvature metric g a b such that g a cgbdP ab p cd /\g\ is spatially constant (as is 
apparent from (7)). For genus 2 and higher the LY equation is much harder to solve. However, we can 
already deduce some properties of Hsd by integrating expression (33), 

H SD = -^(W 2 - 4A ) -R+^(p 2 e- 2 n- (34) 

We have chosen our second class constraints (31) in such a way that the solution /i will not depend on (tt) 
(or A), and therefore the same holds for the last term in (34). Hence, our choice is special in that it produces 
a Hamiltonian that is quadratic in the momentum (tt) conjugate to V. 



3.2 Large Volume Asymptotic Expansion 

Although we can not solve (33) explicitly, our modified LY equation does allow for an interesting perturbative 
expansion. Indeed notice that the volume V appears explicitly in (33) and should be treated as a parameter 
when solving the equation. Therefore we can try to find solutions [i and H$d expanded in powers of 1/V 
and construct the Shape Dynamics in the infinite volume limit. 5 To do this we make the ansatz 

oo oo 

e 2 ^ = Y,n k V- k and H SD = Y^ H k V~ k . (35) 

k=0 fc=-l 

Indeed, from carefully looking at equation (33) it follows that higher powers of V can not occur. From the 
normalization (e 2 '') = 1 we get the restrictions (firj) = 1 and (ilk) = for k > 0. 

The leading order of equation (33) is proportional to V and fixes iJ_i = — \{{^) 2 — 4A). At order V° 
the equation then reads 

- n Q H -R + 2A ln(fio) = 0, (36) 

which is clearly solved by Hq = — R and Hq = 1. The LY equation now becomes 

1 l P 2 -(l + ^ + ...)(-R+^ + ^ + ...)-R + Aln(l + ^ + ---)=0. (37) 



Vl + ?h- + ■■■*' V" ' v ' J V " ' v ' V 2 ' J "V ' V 



If we define the polynomials A k and B k in fii through ft k by 



o 



"Z^^TiT and ln (! + Z. vk) - Z. Vk = v 2- "W ( 38 ) 



1 + V % ^ V k ^ V k ^ V k V *-^ v k 

then from equation (34) or from integrating (37) it follows that for k > 1 

H k = (p 2 A k _ 1 [Q,]). (39) 

The order V~ k equation in (37) then allows us to solve Q k explicitly in terms of Oo through Ofc-i, 



(A + R)- 1 I -A k ^[n}p 2 - ABh [ft] + J2 Hiftk-i ) , 



n^iA + R)- 1 l-^ fc _i[n]p 2 -AB fc _i[n] + 5jj?',nfc_i I, (40) 

where the operator A + R is negative definite and therefore has a well-defined inverse. 

We have therefore obtained a general algorithm to solve the modified LY equation order by order through 
the recurrence relation (40) together with (39). We have calculated the first few H k explicitly leading to a 
Hamiltonian 

Hsd = ~^((n) 2 ~4A)~R+±(p 2 ) + ^((p 2 -(p 2 ))(A + R)-i(p 2 -(p 2 ))) 

+5 ^ [fc(((A + i?)- V - (P 2 ))) 3 ) + 3 {(p 2 + (P 2 »((A + RyHp 2 - (P 2 ))) 2 )} +■■■■ 

In general H k will be a function homogeneous in p 2 of order k. 7 

The expansion has features similar to a tree level (Feynman diagram) expansion with propagator (A + 
R)^ 1 and source term p 2 /V. To make this connection more explicit, let us view the modified LY equation 
(33) as the Euler-Lagrange equation of some action Sly- Such an action can be easily constructed, 

Sly\m, H]=J d 2 x^g L&p R„ - (H - |) (e 2 » l) - \^P 2 ^ , (42) 



5 Notice that this can not easily be done in the reduced phase space approach, in which also the volume itself has to be solved 
for in terms of g and the momenta. 

6 The first few polynomials are A = 1, Ax = -Clx, A 2 = -U 2 + O 2 , A s = -Cl 3 + 2Hin 2 - fif and B = 0,-Bj = -ft?/2, 

b 2 = nf/3-ain 2 . 

7 As the H^ are actually functions on the cotangent bundle to Teichmiillcr space, one might ask how natural they are from 
the perspective of Tcichmuller spaces. As a partial answer we notice that Hi = (p 2 ) is related to the canonical Weil-Pctcrsson 
metric [15, 17], while H 2 is closely related to its curvature [16]. 



where \i is now viewed as an unrestricted function since the Lagrange multiplier H = ^{H$d + ^■(( 7T ) 2 ~ 
4A) + R) enforces the constraint (e 2fi ) = 1 on variation. We can rewrite (42) by singling out the quadratic 
part in fi and H, 

Sly[v,H] = J d 2 xVg L(A + R)fi - 2 t iH + R(pi 3 + ^ 4 + . . .) - ff(2/i 2 + ...)- \^P*) ■ (43) 

The Feynman rules can be read off and we can find H (and therefore H$d) by computing its tree-level one- 
point function. Notice that the action (42) which we use to derive the Hamiltonian is similar to that of two 
dimensional Liouville gravity [18] . More precisely, it is of the form of a Liouville action plus a perturbation 
given by a source term proportional to p 2 /V. 
Two remarks are in order: 

1. We observe from equation (34) that in the limit V —)■ oo the Shape Dynamics Hamiltonian again ap- 
proaches a form that is equivalent to (27). The Hamiltonian is thus asymptotically local and provides 
the volume-changing generator of conformal transformations, so full conformal invariance is asymptot- 
ically attained. 

2. The first three terms in the large volume expansion (41) sum up to an expression equivalent to the 
temporal gauge Hamiltonian S(N = 1), which is local. 

4 Dirac Quantization in Metric Variables 

To expose the difference between Shape Dynamics and General Relativity, we consider the Dirac Quantization 
of pure gravity on the torus in 2+1 dimensions in metric variables, usually referred to as Wheeler-DeWitt 
quantization. For the sake of completeness we first follow [13] and revisit the problems associated with the 
nonlocality arising from the solution of the diffcomorphism constraint in the Whcclcr-DcWitt approach. 
Subsequently we show how these problems are solved by trading local Hamiltonian constraints for local 
conformal constraints, which allows us to perform a Dirac Quantization program for Shape Dynamics. 

4.1 Dirac Quantization of General Relativity on the 2+1 Torus 

Contrary to first order variables, one can not readily quantize General Relativity in metric variables even in 
2+1 dimensions on the torus and sphere. The reason is well explained in [13], which we follow here. Using 
the standard decomposition of the metric and momenta, we can solve the diffcomorphism- constraints for 
the transverse part of the momenta 

^=4( A+ Cf 1(e2A ^ (e ~ 2A7r)) ' (44) 

where k = for the torus. To perform a Wheeler-DeWitt quantization reference [13] chooses a polarization 
for which the configuration operators are given by functionals of the spatial metric and formally considers a 
Schrodinger representation on wave functions ip[f, A;r), which reduces to ip[\;r). Assuming that the inner 

LA. 
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product is constructed from a divergence-free measure, we can quantize the momenta by replacing tt ' 



and p ab — > —i-J?—. The diffcomorphism constraint on the torus still acts nontrivially on the conformal factor, 
and its solution can be quantized as 

This expression is plugged into the Hamiltonian and leads to nonlocal terms in the Wheeler-DeWitt equation 
that are not practically manageable and lead to notorious difficulties in the Wheeler-DeWitt approach [19]. 
To make a connection with reduced phase space quantization reference [13] assumes a solution ?/> [A;t) to 
the Wheeler-DeWitt equation and restricts it to constant York time T through 

^ (T, t) := / V\e l T J d ' x ^ ^ [A; r) (46) 



and inserts this into the Wheeler-DeWitt equation (with vanishing cosmological constant). In terms of the 
Teichmiiller Laplacian A := — rf (c^ + d 2 2 ) this yields 

(( T ^) 2 + A °) MT lT ) = J D\e* T i d2 ^ X (T 2 (e 4A - V 2 ) + 4e 2A AA) ^[A;r), (47) 

where the RHS vanishes if ip has support only on spatially constant conformal factors, while the LHS is 
equivalent to a reduced phase space quantization. 

4.2 Dirac Quantization of Shape Dynamics on the 2+1 Torus 

We now follow essentially the same strategy as in the previous subsection but for Shape Dynamics. Wc 
choose a polarization where functionals of the metric are configuration variables and formally consider a 
Schrodinger representation on functionals ip[X,f;r), such that functionals of the metric are represented by 
multiplication operators and also re-use the representation of the momentum operators given in the previous 

section. We start with the local conformal constraint ir(x) — v U J d 2 yir(y) = 0, which can be readily 

quantized as 



i ( 8 v^O • , j2 



2\5X( X ) V J dy JW)) mf;T ^^ ^ 

where wc work in a chart where the components of g a b are constant. The solution to this constraint is that tp 
depends only on the homogeneous mode of A. We can thus write the general solution to the local conformal 
constraints as a wave function of ip[f; V, r), where V denotes the spatial volume. We now turn to the spatial 
diffcomorphism constraint. Exponentiating the spatial diffeomorphism constraint to finite diffeomorphisms 
implies that for each small diffcomorphism f there is a unitary operator acting as the pull-back under a 
diffeomorphism: 

t// V:[A,/;r)H>V[/ *A,/ o/;r), (49) 

where we assume that the measure is diffcomorphism invariant, so Uf is indeed unitary. The pull-back 
action f*X on the conformal factor is the source of the nonlocal terms that we encountered in the action of 
the diffeomorphisms in the previous subsection. This action is however trivial on the space of solutions to 
the local conformal constraint, since /* V = V. We can thus easily impose the diffcomorphism constraint 

U fo ip = ip (50) 

for all diffeomorphisms f , which implies for solutions to the local conformal constraint that ip[f;V,r) is 
independent of /. We thus find that the solution space to the local constraints of Shape Dynamics consists 
of Schrodinger wave functions ip(V, r). We have not specified the kinematic Hilbert-space that we used, since 
we returned to reduced phase space Hilbcrt space by solving the linear constraints at the quantum level. 8 

We now consider the Shape Dynamics Hamiltonian Hsd = t% (p 2 +P 2 .) — V {( 7r ) 2 ~~ ' 4A), which can be 
quantized on the factor T-L that remains after solving the linear constraints by replacing pi — > —i-J^j and 
(ir) — > —i-M-f. This leads to the quantum Shape Dynamics Hamiltonian 

H=-r 2 (d 2 Tl +d 2 2 )+V 2 (d 2 v + 4A). (51) 

This is the covariant reduced phase space Hamiltonian [10]. We thus confirmed the expectation of the 
last subsection that Dirac Quantization of Shape Dynamics should be equivalent to reduced phase space 
quantization. 



8 We expect that a rigorous construction of the kinematic Hilbert space is possible as follows: First we specify a measure 
dfi on R 3 and consider the Hilbert space "Ho = L 2 (M. 3 ,dfi)/a.e. of square integrable functions modulo the set of functions that 
have support only on sets of measure zero. Next, we construct a Hilbert space W^ as the closure of ffi/gDi//.^/, where Diff. 
denotes the group of diffeomorphisms that are not isometries of reference metrics and where all Hf have one complex dimension. 
Then we construct a Hilbert space T-L c as the closure of ©AgVPCT^; , where A is a conformal factor that preserves the total 
volume of all standard metrics (this is possible since the volume element of the standard metrics does not depend on the point 
in Teichmiiller space of the 2-torus) and where again each ti? has one complex dimension. We define our representation now 
on the tensor product "H = Wo <8> Hd <8> "H c . Informally speaking the Hilbert spaces T-Ld and ~H C consist of the closed span 
of normalized eigenfunctions of the diffcomorphism resp. conformal degrees of freedom of the metric. The reason for this 
construction is that we want to solve the diffcomorphism and local conformal constraints by group averaging, which collapses 
the space "Hd ® H c to C. Solving these constraints thus leads to the Hilbert space M». 
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Excursion: Large Volume Behavior of the Quantum Theory 

We saw in the classical theory that full conformal invariance is attained in the large volume regime. Let us 
now investigate this question in the quantum theory. This is particularly interesting in light of a large-CMC- 
volume/conformal field theory correspondence. For a first investigation we neglect the issue of modular 
invariance and simply investigate the asymptotic volume dependence of solutions ?/>(V, Ti,T2) to the Shape 
Dynamics Wheeler-DeWitt equation 

(-r 2 2 (d 2 T1 + d 2 2 ) + V 2 (d 2 v + 4A)) m n,T 2 ) = 0. (52) 

We can use a separation ansatz ip(V,Ti,T2) = v(V)m(ri,T2) and introduce separation constants a. This 
implies t% (<9^ + d 2 2 ) to(ti,t 2 ) = am(Ti,T 2 ) and 

v "(V) + (4A-^)v(V) = 0, (53) 



which is solved by v±(V) = \JV-K\fke ±il G+v^+V*) (j ^——(2y/AV) ± iY^-—(2VAV)\ in terms of 

Bessel functions J v , Y v . We are interested in 
linearly independent asymptotic solutions are 



Bessel functions J V ,Y V . We are interested in the limit V — ¥ oo, where the term S? vanishes and the two 



«± (V) = e ±i2VXv (l + O (yj } for V -»• oo, (54) 

which is at leading order independent of the separation constant a and hence asymptotically true for all 
solutions to the Shape Dynamics Wheeler-DeWitt equation. The general asymptotic scaling under global 
conformal transformations C = —iVdy is thus 

CVfc(V) = ±2VXVv±(V) fory^oo. (55) 

The difference between the Shape Dynamics wave function and the reduced phase space quantization wave 
function is that the Shape Dynamics wave function is a function of the full metric tp[g] that is constrained 
to be independent of the local degrees of freedom (by diffeomorphism- and local conformal invariance) and 
is thus completely specified by its dependence on V, t\ , r 2 , while the reduced phase space quantization wave 
function is a function i/?(V, ti, T2) on reduced phase space only. We thus have the asymptotic scaling of a 
generic Shape Dynamics wave function under global conformal transformations 

C\j}±\g\ = ±2VXVip±[g] for V -> 00, (56) 

which can be combined with manifest invariance of ip [g] under volume preserving conformal transformations 
Cv PCT{x)'>p±[g\ = to yield the asymptotic scaling of ip[g] under arbitrary conformal transformations 

C(x)^ ± \g]^±2VXy/\g\(x)^ ± [g] iovV^oc. (57) 

Notice the explicit deviations from this scaling at order 1/V, which can be derived by straightforward 
scries expansion of v±(V). Compare this with [20] where a similar result is discussed in the context of the 
AdS/CFT-corrcspondcncc and [21] for a large CMC-volume/CFT discussion. 

5 Possible Consequences for Shape Dynamics in Higher Dimen- 
sions 

Pure gravity in 2+1 dimensions does not possess local degrees of freedom, which is the technical reason for 
many explicit constructions that can not be performed in higher dimensions. However, having an example 
that can be worked out explicitly can give guidance for the treatment of higher dimensions. It is the purpose 
of this section to explore what consequences can be drawn from Shape Dynamics in 2+1 dimensions, in 
particular from the torus universe. A natural generalization of the simplifications occurring on the 2+1 
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torus universe to 3+1 dimensions is the strong gravity or BKL-limit, where spatial derivatives in the ADM 
Hamilton constraints can be neglected, so 

TT ab {x)G abcd {x)TT cd { X ) 



sftfo) = " ^pgy w + 2Av ^y. (58) 

vlslO) 

In this case we can algebraically solve the analogue of the Lichnerowicz-York equation and (e 6( ^) = 1 for the 
Shape Dynamics BKL Hamiltonian in 3+1 dimensions 

Hsd L = V (U^g ac g bd a^) 2 - i(^) 2 g + 2k) , (59) 



where a ab = ■K ab — ^irg ab . We observe that the ingredients V, (y / Tr(a.a)} and (ir) are each invariant 
under volume preserving conformal transformations as well as invariant under spatial diffcomorphisms at the 
expense of nonlocality. The absence of spatial derivatives is what allows us to explicitly construct the BKL- 
Shapc Dynamics Hamiltonian in higher dimensions, i.e. for a similar reason for which the Shape Dynamics 
Hamiltonian could be explicitly constructed for the spherical- and torus- universe in 2+1 dimensions. 9 

In an effort to include spatial derivatives, we could decompose a 3-dimcnsional metric analogously into 
a conformal factor, diffcomorphism and reference metric to the decomposition in 2 dimensions. For this we 
observe that locally one can specify a diffcomorphism class of a metrics by giving three independent curvature 
invariants e.g. (4>i , fa , fa) — [R,R a bR ab , tt ) and integration constants r. Using the Yamabe problem on 

a compact manifold without boundary, one can then impose R(x) = (R) as a gauge fixing condition for the 
conformal factor and thus locally find reference metrics g a b[(R), fa, fa,T; x). Despite this being a purely 
formal construction, since in contrast with the 2-torus and sphere the construction of reference metrics is 
not feasible as it requires the inversion of a complicated system of coupled partial differential equations, one 
finds that one can still not solve for the Shape Dynamics Hamiltonian in all phase space, because of the 
position dependence of g , a , ab (x). Let us therefore consider the restricted phase space 

T,. = f (ff.Tr) G T ADM : ^±(x) = (^°*) and R(x) = (R)\ . (60) 

I 151 \g\ J 

Since we can use R(x) = (R) =: R as a gauge fixing for the conformal gauge symmetry on a compact 
manifold without boundary, we see that the condition a , <J , ab (x) = ( - , <J , ab ) constrains one local physical 
degree of freedom, so F r contains only 3/4 of the physical degrees of freedom of Shape Dynamics. The fact 
that the solution to the modified Lichnerowicz-York equation is homogeneous if all coefficient functions are 
homogeneous and that the volume preservation condition for the conformal factor implies that a homogeneous 



conformal factor is <p{x) = 6 lets us find 



n ab G abcd Tr cd 



H SD \r r = / d 3 x ^ -(R- 2A) ^\{x) , (61) 



i.e. the ADM-Hamiltonian with homogeneous lapse. Observing that on T r we have J d 3 x- — Ga ^ T: 

2 



V ( \J o ah g ac gbdP cd ) ^ T"( 7r ) 2 an d J d 3 Xy/\g\R = VR , where R denotes the Yamabe constant expressed 

as a nonlocal functional obtained by solving R[e ix g;x) = (R) for X [g;x) and R [g;x) := R[e AX ^ 9:x ^g;x), we 
can immediately extend this to the conformal orbit T c of T r to yield 

H SD \ ra = V ((Vv ab 9acg bd <T cd ) 2 g - \{*) 2 - (Ro) + 2A) , (62) 

which differs from the Shape Dynamics Hamiltonian for 1/4 of the physical degrees of freedom. In summary, 
we see that one has to resort to approximation schemes to solve for the classical Shape Dynamics Hamiltonian 



9 The Lichnerowicz-York equation on the torus is algebraically solvable by applying the maximum principle, while the BKL 
Hamiltonian is algebraic to begin with. 
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that is equivalent to General Relativity, which is an obstruction to the direct construction of a quantum 
theory of Shape Dynamics. 

One may understand this as a hint to use the BKL-Hamiltonian as a starting point for Shape Dynamics 
and introduce spatial derivatives as perturbations in the classical theory. Similarly, one may hope to quantize 
the ultralocal theory and develop a perturbation scheme at the quantum level. The feasibility of this idea 
remains to be explored, but if possible it would lead to an expansion of the quantum theory that is exact 
precisely where quantum effects are expected to be most important. 10 

A different direction could be to follow the original construction of Hamiltonians in Loop Quantum 
Gravity, where one first constructs a rigging map to solve the diffcomorphism constraint and subsequently 
constructs a quantization of the General Relativity Hamiltonian constraints on the solutions to the diffco- 
morphism constraint to avoid problems with anomalies. An analogous strategy for Shape Dynamics is to first 
construct a rigging map to solve the local conformal constraint and subsequently construct a quantization 
of the Shape Dynamics Hamiltonian on the solution space of the local conformal constraints. 

6 Conclusions 

The true value of pure gravity on a torus in 2+1 dimensions is that it is a nontrivial yet completely solvable 
model that exhibits many of the features of more complicated gravitational systems. It has thus established 
itself as a valuable testing ground for new gravitational theories such as Shape Dynamics that one can use 
to learn about the new theory. The main difficulty in constructing Shape Dyamics is to obtain explicit 
expressions for the Shape Dynamics Hamiltonian on the full ADM phase space, so we are mainly interested 
in obtaining an explicit Shape Dynamics Hamiltonian. In these investigations we observed 

1. The explicit (rather than formal) solvability of the initial value problem on the 2-torus (and 2-sphere) in 
CMC gauge is the technical reason for the explicit constructability of the Shape Dynamics Hamiltonian 
on these topologies. We find that the Shape Dynamics Hamiltonian formally coincides with the reduced 
phase space Hamiltonian and that this Hamiltonian is invariant under diffeomorphisms and conformal 
transformations that preserve the total volume. The difference between the two is however that the 
Shape Dynamics Hamiltonian is a function of the full ADM phase space that happens to functionally 
depend only on the image of reduced phase space under the canonical embedding, while the reduced 
phase space Hamiltonian is a phase space function on reduced phase space itself. 

2. Although we cannot explicitly construct the Shape Dynamics Hamiltonian on higher genus Riemann 
surfaces, we can construct it perturbativcly. In particular we use an expansion that becomes good in a 
large volume regime. The leading orders of this expansion then turn out to coincide with the temporal 
gauge Hamiltonian. 

3. The Hamiltonian is in general a nonlocal phase space function. However, it becomes a local phase 
space function in the large volume limit; in particular one finds that the leading order in a large 
volume expansion turns the Hamiltonian into the conformal constraint that changes the total volume, 
so full conformal invariance is attained in this limit. 

4. Since all local constraints are linear in momenta, one can quantize these as vector fields on configuration 
space. Then gauge invariance implies that the wave function is invariant under the flow generated by 
these vector fields, which in turn implies that the wave function depends only on reduced configuration 
space, which is finite dimensional. The Hamiltonian depends only on operators that preserve the 
reduced phase space, and thus Dirac quantization of the field theory is effectively reduced to reduced 
phase space quantization. 

5. The construction of Shape Dynamics on the torus relied on the explicit solvability of the modified 
Lichnerowicz-York equation. This suggests to take the strong gravity BKL limit, where the analogous 
equation can be explicitly solved as a starting point for quantization, and to introduce the effect of 
derivative operators through a perturbation series. However it may turn out that alternative quanti- 
zation strategies are more advantageous as, e.g., the one discussed at the end of section 5. 



10 The BKL-conjecture [22, 23, 24] roughly states that spatial derivatives are negligible when a generic singularity is ap- 
proached, which is the regime where quantum effects of gravity are expected to be most important. 
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Lastly, let us remark that Shape Dynamics is a natural setting to discuss a quantum "large-CMC- volume/CFT" 
correspondence, which we were able to investigate explicitly in an excursion at the end of section 4, where 
we found an explicit asymptotic scaling of solutions to the Wheeler-DeWitt equation under conformal trans- 
formations, which is very similar to the correspondence explored in [20]. 
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